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Abstract 

A simple pseudopotential model is proposed, which allows the phonon spectra and tem- 
perature dependence of the lattice heat capacity of Ir and Rh be described with a high enough 
accuracy. A careful comparison of the calculated and experimental values of the lattice heat 
capacity is carried out, with the procedure of the identification of the phonon contribution to 
the heat capacity and determination of the characteristics (momenta) of the phonon density of 
states from the experimental values of the total heat capacity of metal at a constant pressure 
being described in detail. The results of the theoretical calculations explain, in particular, 
such peculiar feature of Ir and Rh, unusual for cubic metals, as a sharp (more than by a factor 
of 1.5) decrease in the effective Debye temperature with increasing termperature. The tem- 
perature dependence of the mean square amplitude of atomic displacements in Ir and Rh has 
been calculated. Basing on the band calculations the manifestation of the Kohn singularities 
in the phonon spectra of Ir are discussed. 

§1. Introduction 

From the physical point of view iridium and its analogue rhodium stand out 
by their unusual mechanical properties among other FCC metals (see discussion 
in Gornostyrev et al 1994). The peculiarities of mechanical properties (at least 
those inherent to high pure single crystals) should be eventually determined by 
the specific features of interatomic interactions in metal. Greenberg et al 1990 
and Ivanov et al 1994 succeded in constructing of a pseudopotential model of pair 
forces for these metals wich describes surprisingly well (for transition metals) 
their elastic properties. In terms of this model the peculiarities of the Ir and Rh 
defect structure (Ivanov et al 1994, Gornostyrev et al 1994) could be analyzed, 
and the qualitative explanation of some of their peculiar mechanical properties, 
primarily, the possibility of brittle failure of single crystals after a long-term 
stage of plastic deformation could be given. At the same time it was found 
that in the description of the dispersion curves of phonons in Ir measured by 
Ivanov et al 1994 the model does not prove to be accurate enough for oscillation 
modes near the Brillouin zone boundaries. According to the calculations of 
Ivanov et al 1994 the maximum frequency of the phonon spectrum appears to 
be overestimated by 25 % comparing with the experimental value found from 
the results of tunnel experiments (Kraynukov et al 1989), and flattenings in the 
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(110), (100) directions are observed, which are not described by the proposed 
model. 

The main purpose of the present paper is to ascertain the physical nature of 
these descrepancies, improve the model correspondingly, and calculate a number 
of thermodynamical properties with the allowance for the corrected phonon 
spectra. 

It may seem that the use of a simple pseudopotential model for these pur- 
poses does not correspond to the current state of the quantitative theory of 
metals, when completely ab initio methods for calculation of phonon spectra 
have been developed (see, e. g., review Maksimov and Savrasov 1995). How- 
ever, the description of a wide range of lattice properites, including temperature 
dependence of different values is still beyond the capabilities of completely ab 
initio approaches. Moreover if it is taken into account that the pseudopotential 
approach also gives the pairwise potentials which can be used for the simulation 
of defects, diffusion processes, radiation effects, for the description of the liq- 
uid phase properties (for example, the maximum value of the structural factor, 
determined by the packing parameter) etc. then the importance of the develop- 
ment of the pseudopotential theory for the metals for which it is valid appears 
indisputable. An important advantage of this approach is also its simplicity 
enabling many complicated problems to be discussed not only at the qualitative 
but also at the quantitative level, and the close relationship with the chemical 
bond language used in the material science (Gornostyrev et al 1994). For the 
clarification of the physical mechanism of the Ir dispersion curve flattening the 
results of band calculations are given in §2, on which basis the Kohn singularities 
in the phonon spectra are discussed. In §3 a modification of the pseudopotential 
form permitting these effects to be partially simulated for the purpose of an op- 
timal description of the phonon spectra is proposed. In §4 the characteristics of 
the phonon spectra extracted from the experimental data on the Ir and Rh heat 
capacity are given with a detailed description of the corresponding treatment 
procedure. §5 presents the calculation results of the lattice heat capacity, and 
their comparison with the experiment. §6 summarizes the results of the paper. 

§2. Fermi surface and Kohn singularities in Ir. 

Generally, one could think that the high accuracy of the description of clastic 
moduli, and hence, of phonon spectra at small wave vectors q, provide a good 
agreement of the theory with experiment for the dispersion curves in the larger 
part of the Brillouin zone. Significant deviations from the "smooth" behavior 
of the phonon frequencies depending on the wave vector q can be connected, 
first of all, with the Kohn singularities (Izymov and Chcrnoplckov 1995) when 
q coincides with the Fermi surface diameter 2kp in the corresponding direction. 
However in the three dimensional case the typical Kohn singularities in the 
phonon frequencies w(q) are of the order (q — 2fcp) In \q — 2fcp| and are weakly 
expressed. The Kohn singularities clearly expressed in w(q) can be only observed 
in the case of the Fermi surface flattening in one (cylindrical sections) or two 
(flat sections) directions when Su)(q) <~ ±\/\q — 2/c_f|9(±((/ — 2fcp)) (6(x > 0) = 
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1, Q(x < 0) = 0) or 5ui(q) ~ In \q — 2kp\, respectively (Katsnelson, Naumov and 
Trefilov 1994). Figs. 1,2 present the results of our calculations of the Fermi 
surface of Ir and nonuniform susceptibility 



where e„(k) is the electron energy spectrum; n, n' are the band indices; /(e) is 
the Fermi function; integration over k goes over the Brillouin zone (calculation 
method see in review by Katsnelson et al 1994). The function x(q) determines 
the form and magnitude of the Kohn singularities in w(q). In fig. 1 the flat- 
tened sections of the Fermi surface in the directions (100) and (110) are clearly 
seen. Fig. 2 shows the maxima of x(q) a t the relevant values of wave vectors 
(in the case of the spherical Fermi surface the maxima should have manifested 
themselves only in dx/dq). At the same time in the calculations of the phonon 
spectra in terms of the pseudopotential model in the second order of the per- 
turbation theory (Greenberg et al 1990, Ivanov et al 1994) the character of the 
Kohn singularities corresponds to the assumption about the sphericity of the 
Fermi surface and, therefore, the phonon spectrum in the relevant region of the 
q space, naturally, proves to be too rigid. If we want to improve the description 
of phonon spectra in terms of the simple pseudopotential approach, the soften- 
ing of the phonon frequencies due to the non-sphericity of the Fermi surface can 
be simulated and described by means of changing the form of the Fourier image 
of pseudopotential V ps (q) at Q.6(2n/a)^q^(2Tv/a) (a is the lattice constant), de- 
termined from the band calculations results. It is this circumstance that should 
be primarily taken into account in changing the pseudopotential model for Ir 
and Rh. 



In the calculation of phonon spectra, elastic moduli, and the interatomic 
interaction potential in papers by Greenberg et al 1990, Ivanov et al 1994 one 
of the common pseudopotential models (Heine and Wcaire 1970; Browman and 
Kagan 1974; Vaks and Trefilov 1977) was used. The Fourier component of the 
local pseudopotential of the electron-ion interaction was taken in the form of 
the Animalu-Heine pseudopotential 



where fio is the atomic volume, Z is the effective ion charge, rg, Vq are the fitting 
parameters determined from the condition of the agreement with experiment 
Qq and shear modulus -B44 at pressure p = 0; £ is the constant whose value 
is selected proceeding from the consideration of fast enough convergence of the 
sums over the reciprocal lattice vectors (here and below the atomic units are 
used). The screening function of the electron gas was taken in the approximation 
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§3. Calculation of Ir and Rh phonon spectra 
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by Geldart and Taylor 1970. See Greenberg et al 1990, Ivanov et al 1994 for 
the details of the theoretical model, including the substantiation of the choice of 
values Z = 4.5 for Ir and Z = 3.86 for Rh. As discussed in the introduction this 
model provides a high enough accuracy of the description of elastic moduli and 
much lower one in the description of the phonon frequencies near the Brillouin 
zone boundaries. The latter is the result of neglecting the nonsphcricity of 
the Fermi surface, which gives no possibility to describe adequately the Kohn 
singularities in the phonon spectra (see § 2). 

For the purpose of simulating these effects remaining in the framework of 
the simplest pseudopotcntial model, we considered various possible changes in 
the form of the psedopotential form factor 

V ps (q)=V p A s H (q) + SV ps (q) (3) 

at large q. It is seen from the results of §2 that the "dangerous" region of the 
wave vectors, containing the strongest Kohn singularities (which is the result 
of the maximum non-sphericity of the Fermi surface in the relevant directions) 
is approximately (0.5 < g\ < 0.7)gi, where g\ is the minimum vector of the 
reciprocal lattice. The addition to the pseudopotential 5V ps (q) to simulate these 
singularities should be noticablc in the relevant region of q values and vanishingly 
small at q = and q > g\. It is natural to expect that the addition SV ps (q) 
would not noticably affect the calculation values of pressure and elastic moduli 
which are solely determined by the behavior of V ps {q) at q — > and by V^g) 
for the reciprocal lattice vectors. The concrete form of the addition 8V ps (q) is 
rather arbitrary, however it should affect the V ps (q) in the region of q specified 
at the end of § 2. The following dependences satisfying the above mentioned 



condition were actually studied: 

SiV ps (q) = (3(qr ) cxp[-a(or ) 4 ] (4a) 

52V ps {q) = p{qr f exp[- a (qr ) 4 ] (46) 

S 4 V ps (q) = /3(qr ) 2 cxp[-a(q - q a frl] (4c) 

S 5 V ps (q) = P(qr ) 2 cxp[-a(q - q ) 4 r 4 ] (4d) 



where parameters /?, a (as well as r , Vb) are found from conditions f^o = ^q 0P , 
Qlt — ®lt' w x — ^tnax- However, exact fitting to these conditions turned 
out to be impossible, so the set of parameters giving the best approximation 
was used. Quantity ujx is the frequency of phonons at the X point of the 
Brillouin zone (this is the point where the phonon spectrum calculated using the 
Animalu- Heine pseudopetential used as a basis for the fitting procedure reaches 
its maximum) . The experimental values of uj max were taken from results of the 
tunnel experiments by Kraynukov et al 1989 for Ir and Zhalko-Titarenko et al 
1989 for Rh. The experimental of Lo ma x and calculated values of ujx are listed 
if table 2 for Ir and table 3 for Rh. q is the wave vector at the point X of the 
Brillouin zone. Among the additions (4a)-(4d) (4c) gives the best agreement of 
the calculated phonon spectra with the experiment, although the choice of the 
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addition if the form (4a) permits nearly the same agreement between the theory 
and the experiment to be reached. It is important to emphasize that parameter f3 
is choosen so that additions (4) were small without giving significant distorsion 
of the V ps (q) form when comparing with the standart Animalu- Heine model. 
Nevertheless, calculation results for the phonon spectrum (but not for q — > 0, 
i. e. not for elastic moduli) proves to be very sensitive to the explicit form 
of these small additions. The sensitivity of the calculation results to the form 
of the pseudopotential will be considered in detail in a separate paper. Table 
1 lists the parameters for two versions of the pseudopotential, V ps H (q) and 
V^ H (q) + SV ps (q), where addition 5V ps (q) is choosen as (4c) which, comparing 
with the rest possibilities investigated gives the best description of the lattice 
properties of Ir and Rh. Tables 2 and 3 and figs. 3,4 give the results of the 
corresponding calculations. Qlt — Bd(T = 0) is the low temparature value 
of the Debye temperature calculated from the theoretical values of Bik, using 
as the bulk modulus B its static value B st determined by differentiating the 
total energy over volume. Quantity Bd yn obtained from the phonon spectra by 
extrapolation has a methodical meaning: difference \B st — Bd yn \ characterizes 
the value of multiion interactions (Browman and Kagan 1974). For Ir and Rh its 
value is substantially smaller than for other transition metals (Vaks et al 1989). 
Quantity 70 = —d In Qlt/9 In fl is the low temperature Gruneisen parameter. 
Note the high sensitivity of 70 to the form of the pseudopotential (see tables 2- 
3). We will discuss the sensitivity of results to the form of the pseudopotential in 
a separate paper. Unfortunately, at present we don't have experimental results 
to compare with the calculated values of 70. The arrows in figs. 3,4 show 
the positions of the Kohn singularities in the approximation of free electrons 
|q + g| = 2/cf- An unusual feature of the considered phonon spectra is that in 
one of the branches in the X-W directions ui(q) is practically independent on 
q. This feature is not connected with the Kohn singularities in any way and is 
obtained even in the approximation of the nearest neighbours for interatomic 
interactions. As discussed in papers by Ivanov et al 1994, Gornostyrev et al 
1994 this approximation describes reasonably the clastic moduli of Ir. 

It is seen that by introducing "simulation" corrections into V ps (q) one can 
reach a better description of the phonon spectra of Ir and apparently Rh. This 
allows us to expect that the calculated phonon spectra can be used for obtaining 
true enough information on the thermodynamical properties of these metals. It 
should be emphasized that it is difficult to obtain the information about the 
phonon spectra of Ir using the method of inelastic neutron scattering because 
of the large cross section of neutron absorption. Therefore when considering 
the agreement between the theory and experiment, it is important to get the 
information on the Ir lattice dynamics from the analysis of the temparature 
dependences of the thermodynamical properties of Ir, primarily its heat capac- 
ity. In the next section we shall use this procedure for obtaining the values of 
different characteristics (momenta) of the Ir and Rh phonon spectra. 
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§4. Analysis of experimental data on Ir and Rh heat capacity 



The lattice heat capacity is determined by the form of phonon density of 
states g(u>). Some characteristics of the latter, which can be used for the com- 
parison of our description of phonon spectra with the experiment, can be derived 
from the analysis of the temperature dependence of the heat capacity at a con- 
stant pressure C P (T). 

In this section we identified the lattice component in the heat capacity mea- 
sured in the experiment and determined values of some average frequencies 
(momenta) of the phonon spectra. 

We used the following set of data in the analysis: in region 12-250 K - 
direct experimental data (Clusius and Losa 1959) and in region 1.2-12 K — the 
data from the handbook by Hultgren et al 1973, using the data by Wolcott 1956 
for the region 1.2-20 K. 

In the analysis we neglected the temperature dependence of the coefficient 
of the electron heat capacity and assumed that the anharmonicity is weak. 
The assumption on the smallness of the anharmonicity contributions to the 
thermodynamical properties of metals is justified as it was verified by direct 
calculations and proved to be true even for such metals with soft phonon modes 
as alkali ones (Vaks, Kravchuk and Trefilov 1980). With these assumptions the 
experimentally measured heat capacity at constant pressure C p is described by 
the relationships 

C p = C ph + C a (5a) 

C a = { 1 +(A- 1 )(C ph /2,R) 2 }T (56) 

where C p is the molar heat capacity at constant pressure; 

C p h is the phonon component of the heat capacity in the harmonic approx- 
imation; 

C a — contains the contributions to the heat capacity, having the linear 
temperature dependence caused by anharmonic effects, thermal expansion of 
the lattice as well as by the conduction electrons; and here 

7 is the coefficient of the electron heat capacity at low temperatures; 

A is the coefficient at the linear temparature term at high termperatures; 

R is the gas constant. 

Interpolation formula (5b) for C a gives the relevant linear temperature asymp- 
totes both at low and high temperatures and ensures a smooth transition be- 
tween the low- and high temperature asymptotes by a law similar to the Nernst- 
Lindeman formula (Rcznitskii 1981). 

In the high temparature region the phonon component of the heat capacity 
was described by the expression proposed by Naumov (Naumov 1994): 



- " eXP{Z) 2 - f * - V + ^ ) (7) 
(l-exp(z)f V 12 240 1 
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Here the asymptotic expansion of the phonon heat capacity by small parameter 
z = fl/T is used: 

where B n are the Bernulli numbers (B2 = 1/6, B4 = —1/30, Bq = 1/42, 
B s = —1/30, B w = 5/66 and for all odd n all B n = beginning from n = 3). 

In (6) corrections of second and fourth order by fl/T are written separetaly, 
and function <p(^r) allows for all higher order corrections in the "Einstein" 
approximation, i. e. setting Ct n — f2» for all n > 6. 

Quantities 0„ characterize the momenta (average frequencies) of the phonon 
spectrum according to the relationship 

roc I poo 

(ftj" =< w ™ >= J g(u;)uj n du J g{uj)duj, (9) 

g{u>) being the phonon density of states. 

Quantity 7 in C a , (5b) was determined by the standart procedure: by the 
approximation of the heat capacity in the low temparature region (particulary 
1.2-15 K) by a relationship such as C p = 7T + f3T 3 + aT 5 . Quantities 7, (3 and 
a as well as the value of the Debye temperature in the low temperature limit 
9lt, related to /3 by the equality /3 — 12tt 4 R/(5Q^ t ) are listed in the tabic 4. 
The estimates of 7 and Olt obtained from the set of all used data (Clusius and 
Losa 1959, Wolcott 1956) agree resonably well with the data available from the 
literature (see table 4). 

As for quantities A, f^, ^4 and f2* they were determined by the method of 
least squares, approximating the heat capacity by relations (6,7) in the tem- 
perature region 50-250 K. Within this temperature range these relationships 
described the experimental results of Clusius and Losa 1959 with a mean-square 
deviation of about 0.6 %. The values of parameters A, SI2, ^4 and f2», deter- 
mined by the least square method are listed in table 4. The table also gives 
the value of the Debye temperature in high temperature limit &ht = ®d(T), 
where @lt < T < T m (T m is the melting temperature) related to the second 
momentum of the phonon spectrum by: £1 2 = ®HTy/3/5 (Maradudin, Montroll 
and Weiss 1963). 

The described analisys permitted phonon contribution C p h in the harmonic 
approximation to be separated from the total heat capacity in the temperature 
region below 273 K, excluding the electron and anharmonic contributions. Some 
momenta of the phonon spectra are expressed directly by the integrals of the 
phonon heat capacity (Junod 1980). 
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< 10 




(10) 



We have calculated these momenta. In the region 15-273 K the integration 
was carried out by the experimental points (Clusius and Losa 1959) while outside 
this region we extrapolated the heat capacities by a low temperature asymptote 
C p h = ftT 3 + aT 5 and high temperature asymptote in the model of the Debye 
spectrum as in Mirmelstein et al 1984. The average frequencies corresponding 
to this momenta are listed in table 4. The meaning of the quantities 0_2, 
f2i corresponds to the meaning of fl n in relationship (9) for n = —2, —1 and 1, 
while the meaning of Qi og is determined by 



§5. Peculiarities of the temperature dependence of heat capacity 



Although the directly measured quantity is heat capacity, it is more con- 
vinient to discuss its temperature dependence via the temperature dependence 
of the effective Debye temperature &(T). This is related to the circumstance 
that the C p h (T) reaches the asymptotic 3R rather quickly with the increase of T 
(and this result will be reproduced in any, even apparently incorrect model of the 
phonon spectrum, and the dependence O(T) is a finer characteristic determin- 
ing, in particular, the rate at which C p h(T) reaches the asymptotic (Maradudin, 
et al 1963). 

An unusual feature of Ir and Rh, comparing with the other FCC metals, is 
the sharp fall of the effective Debye temperature with temperature (0lt/0ht~1-5). 
For example according to handbook by Zinoviev 1989 for Ca 9^t ~ ®ht ~ 
230 K, for Cu 9 LT « 342 K, Q HT w 310 K, for Pd 6 LT « 283 K, Q HT w 
275 K. Only in Ni ratio & lt /®ht « 476 if/345 K w 1.4 differs strongly enough 
from unity (note, by the way, that in this handbook Ir has a very underestimated 
value of Qht)- 

For comparing the experimental data on C p h (T) in Ir and Rh the calculations 
were made using phonon spectra w q „ that we had obtained 



where the summation over the wave vector is made over the Brillouin zone and 
v is the number of the branch. The effective Debye temparature was derived 
from the standart equation (Maradudin, et al 1963). 




of Ir and Rh 




(11) 




0(T)/T 



(12) 
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The calculation results are shown in figs. 5,6 and in table 4. It should be noted 
that the corrected pseudopotential with the addition of the form (4c) describes 
average frequencies 0„ much better than the Animalu-Hcinc pseudopotential 
proposed by Ivanov et al 1994. For example, the agreement between the the- 
oretical and experimental values of fl± for Ir improved from 17 % to 9% and 
similarly for the rest Q n . It should be emphasized that in the calculation of 
CV(T) the change in the volume with temperature was not taken into account. 
It may be assumed that taking into account this dependence could improve 
the agreement between the theory and experiment by some more percent. The 
difference between the experimental and theoretical values of effective Debye 
temperature (\(Q e xp — @theor) /®theor\) does not exceed 10 %. It seems sur- 
prising that a simple pseudopotential model can provide the description of the 
lattice heat capacity of transition metals (Ir and Rh) at an accuracy level com- 
parable with such for alkali metals. The fact is really unusual as it is commonly 
believed that the theory of pseudopotential in its simplest form used in this pa- 
per is a priori not applicable to transition metals and has not so high accuracy 
even for polyvalent s — p metals (e. g. Al). Fig 7 shows the corresponding data 
on 6_d(T) for Na obtained from the calculation results (Vaks et al 1978). The 
accuracy of the description of the phonon spectra in the pseudopotential model 
for other transition metals is incomparably lower while, say, for Ni this model 
is totally unsuitable (Vals et al 1989). 

When phonon spectrum aj q „ is known in the whole Brillouin zone one can 
calculate the temperature dependence of the average square of the atomic dis- 
placements in the harmonic approximation 

<13) 

where M is the atom mass. The corresponding results are shown in fig. 8. As 
far as we know the experimental data on x 2 {T) in Ir and Rh are currently not 
available. 

§6. Conclusion 

In conclusion we shall underline the most important results obtained in the 
paper. The most surprising and interesting finding appears to be that the 
simplest physical pseudopotential model, we proposed earlier, with appropriate 
modifications seems to provide a high accuracy of the description of phonon 
spectra and lattice heat capacity of transition metals Ir and Rh and the accu- 
racy of description is comparable to that for alkali metals. This confirms our 
earlier consideration (Greenberg et al 1990, Gornostyrev et al 1994) about the 
applicability of the model of pair interactions for the description of the lattice 
properties of Ir and Rh. However the effects of non-sphericity of the Fermi sur- 
face, substantial for all transition metals cannot be completely accounted for. 
Nevertheless in Ir (and, as appears in Rh) the role of these effects is reduced 
to some distorsions of the despersion curves in directions in the Brillouin zone 
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((100) and (110)) and appears to be insignificant in the integral characteristics, 
i. e. in the heat capacity. This situation would be absolutely impossible if in 
these metals the Fermi level were at the peak of the density of states Ep, as, 
say in Ni or Pd, because this peak affects substantially the elastic moduli and 
phonon frequencies in the noticable part of the Brillouin zone (Katsnelson et al 
1994). The band calculations show, however (Greenberg et al 1990) that not 
only peak of N(E) does not exist but also there are no relatively weak van Hove 
singularities near Ep in Ir. Of course in itself it cannot be taken as a direct 
evidence of the applicability of the psudopotential model (in W, say, N(E) near 
Ep is also a smooth function but a reasonable model pscudopotential cannot 
be constructed). However this can be considered as a some indication to the 
possibility of such simple description. In the whole, to solve the question of its 
applicability a wide range of the lattice properties (thermal expansion, equation 
of state at finite temperatures, other anharmonic effects, etc.) has to be studied. 
We shall consider these questions in separate papers. 

The study described in this paper became possible partially due to the fi- 
nancial support of the International Science Foundation (Grant RGQ300) and 
Russian Fund for Fundamental Research (Grant 95-02-06426). 
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Captions to tables of the paper by Katsnelson et al "Peculiarities 

OF PHONON SPECTRA ..." 



Table 1. Pseudopotential parameters (in atomic units) in relations (2), (4). 

Table 2. Values of energy E, elastic moduli Bik, their derivatives over pres- 
sure, Debye temperature at T — > Qlt, Gruneisen parameter 70 and maximum 
phonon frequency uj max for Ir. E — Ry/at, P,B — 10 11 dyn/cm 2 , 6lt — K, 
Umax — THz. AH corresponds to the values calaulated by the Animalu-Heine 
pseudopotential, 4c — the same with the addition (4c). 

Table 3. Values of energy E, elastic moduli their derivatives over pres- 
sure, Debye temperature at T — > (@lt), Gruneisen parameter 70 and maxi- 
mum phonon frequency for Rh co max • E — Ry/at, P,B — 10 11 dyn/cm 2 , @ LT 

— K, UJ m ax THz. 

Table 4a. Parameters characterizing the experimental heat capacity at a 
constant pressure, low temperature and high temperature limits of the Debye 
temperature and the phonon spectrum momenta for Ir. AH corresponds to 
the values calculated by the Animalu-Heine pseudopotential, 4c — the same 
with the addition (4c). [Clu] and [Geb] are data sets from Clusius et al 1959 
and Geballe et al 1966 respectively Subcolumn "this paper" in the column 
"experiment" containes values derived from the experimental data using the 
technique described in § 4. 

Table 4b. Parameters characterizing the experimental heat capacity at a 
constant pressure, low temperature and high temperature limits of the Debye 
temperature and the phonon spectrum momenta for Rh. AH corresponds to the 
values calculated by the Animalu-Heine pseudopotential, 4c — the same with 
the addition (4c). [Clu], [Geb] and [Wol] are data sets from Clusius et al 1959, 
Geballe et al 1966 and Wolcott 1956 respectively. Subcolumn "this paper" in 
the column "experiment" containes values derived from the experimental data 
using the technique described in § 4. 
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Figure captions to the paper by Katsnelson et al "Peculiarities of 

PHONON SPECTRA ..." 



Fig.l Section of the Fermi surface of Ir. 

Fig. 2 Non-uniform susceptability x(q) m directions (110) (a) and (110) (b). 
Numbers in fig. 2b designate contributions from the corresponding branches. 

Fig. 3 Phonon spectrum in Ir. Dashed line — the phonon spectrum calculated 
by the Animalu-Heine pseudopotential, solid line — the same with the addition 
(4c). Experimental points (o) are taken from Ivanov et al 1994. Arrows indicate 
the Kohn singularities in the approximation of free electrons. 

Fig. 4 Phonon spectrum in Rh. Dashed line — the phonon spectrum calcu- 
lated by the Animalu-Hcinc pseudopotential, solid line — the same with the 
addition (4c). Arrows indicate the Kohn singularities in the approximation of 
free electrons. 

Fig. 5 Dependence of the lattice heat capacity at a constant volume (a) and 
the effective Debye temperature (b) on temperature for Ir. Dashed line corre- 
sponds to the values calculated by the Animalu-Heine pseudopotential, Solid 
line — the same with the addition (4c) . The experimental values derived in this 
paper denoted as o. T p i = huj p i = h(4:TrZe 2 /Mflo) 1 ^ 2 , u p i is the ionic plasma 
frequency. 

Fig. 6 Dependence of the lattice heat capacity at a constant volume (a) and 
the effective Debye temperature (b) on temperature for Rh. Dashed line cor- 
responds to the values calculated by the Animalu-Heine pseudopotential, Solid 
line — the same with the addition (4c) . The experimental values derived in this 
paper denoted as o. T pi = fno p i = H(4:TrZe 2 /MQq) 1 / 2 , uj p i is the ionic plasma 
frequency. 

Fig. 7 Theoretical ( — ) and experimental (o) temperature dependences of the 
lattice heat capacity at a constant volume (a) and the effective Debye temper- 
ature (b) for Na from Vaks et al 1978. T p i = hui p i = h{AirZe 2 /Mft ) 1/2 , uj p i is 
the ionic plasma frequency. 

Fig. 8 Average square of atomic displacement in Ir (a) and Rh (b). Dashed 
line corresponds to the values calculated by the Animalu-Heine pseudopotential, 
Solid line — the same with the addition (4c). T pi = hui p i = Ti(<iirZe 2 /MVLq) 1/2 , 
u p i is the ionic plasma frequency. 
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table 1 



metal 


n 


Z 




ro 


V 


a 


(3 


Ir AH 


95.88 


4.5 


0.3 


2.696 


-1.243 






Ir 4c 


95.88 


4.5 


0.3 


2.696 


-1.243 


1.285 


-4.29-10 -3 


Rh AH 


92.57 


3.86 


0.3 


2.678 


-1.251 






Rh 4c 


92.57 


3.86 


0.3 


2.678 


-1.251 


1.276 


-6.76-10" 3 



table 2 





E 


P 


B s t 


Bdyn 


#33 




8P 


dB 44 
BP 


dB 3t 
8P 


dBdyn 

8P 


Qd 


7o 




AH 


-9.064 


0.00 


31.27 


35.17 


14.01 


26.69 


2.50 


4.63 


5.95 


5.90 


409 


1.73 


8.06 


4c 


-9.060 


0.00 


29.26 


33.14 


14.40 


24.43 


2.21 


5.77 


6.71 


6.60 


401 


1.99 


6.98 


exper 




0.00 


37.1 


37.1 


17.1 


26.9 










425 




6.91 



table 3 





E 


P 


B st 


-E^dyn 


-B33 


-B44 


dB 33 

dp 


<9B 4 4 
3P 


dB Bt 
dP 


t>B dyn 
3P 


e D 


7o 


^rnaxiX) 


AH 


-6.860 


0.00 


25.33 


29.02 


10.97 


22.15 


2.09 


4.30 


5.70 


5.51 


500 


1.54 


10.03 


4c 


-6.857 


0.00 


22.30 


25.91 


11.49 


18.76 


1.60 


6.40 


7.21 


6.87 


482 


2.05 


7.70 


exper 




0.00 


26.9 


26.9 


11.5 


19.4 










482 




7.63 
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table 4a 



parameter 


experiment 


theory, this paper 




mini 


1 V-TCU 1 


AH 




7, mJ mol" 


l K -2 


3.26 


3.18 


3.20 






(3, mJ mol" 


-l K -4 


0.0228 


0.0244 








a, mJ mol" 


-l K -6 


0.00005 










e LT , k 




440 


430 




409 


401 


Qht, K 




290 


289 




355 


321 


A, mJ mol" 


-l K -2 


4.7 














192 






227 


214 


0_ 2 , K 




188 






222 


213 


n-i, k 




202 






240 


224 


fii, K 




217 






260 


241 


n 2) k 




225 






267 


248 


4 , K 




236 






276 


261 


a,, k 




240 











table 4b 



parameter 


experiment 


theory, this paper 


this paper 


[CluJ 


[Budj 


[Wolj 


AH 


4c 


7, mJ mol" 


l K -2 


4.95 


4.18 


4.65 


4.90 






0, mJ mol" 


-l K -4 


0.0195 


0.0213 










a, mJ mol" 


-l K -6 


0.00002 












e LT , k 




505 


450 


512 




500 


482 


Qht, K 




345 


346 






437 


362 


A, mJ mol" 


-l K -2 


5.5 












fllog, K 




227 








283 


246 


0_ 2 , K 




222 








281 


245 


K 




240 








299 


256 


fti, K 




260 








328 


273 


n 2 , k 




267 








339 


280 


4 , K 




276 








359 


293 


a., k 




278 
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